Local interactions and non-Abelian quantum loop gases 
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Two-dimensional quantum loop gases are elementary examples of topological ground states with Abelian or 
non-Abelian anyonic excitations. While Abelian loop gases appear as ground states of local, gapped Hamil- 
tonians such as the toric code, we show that gapped non-Abelian loop gases require non-local interactions (or 
non-trivial inner products). Perturbing a local, gapless Hamiltonian with an anticipated "non-Abelian" ground- 
state wavefunction immediately drives the system into the Abelian phase, as can be seen by measuring the 
Hausdorff dimension of loops. Local quantum critical behavior is found in a loop gas in which all equal-time 
correlations of local operators decay exponentially. 
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Non-Abelian topological phases are the focus of considerable 
excitement as a result of their universality, their novelty, their 
beautiful mathematical properties, and their potential applica- 
tion to quantum computing iH. However, the only concrete 
physical system in which there is any experimental evidence 
for a topological phase is the two-dimensional electron gas at 
high magnetic fields, i.e. in the quantum Hall regime. In order 
to find such phases elsewhere - in transition metal oxides or in 
ultra-cold atomic gases - it is important for theory to serve as 
a guide by identifying conditions which a system must satisfy 
in order to support a non-Abelian topological phase. 

One simple class of models is associated with quantum loop 
gases (QLG), in which an orthonormal basis of the low-energy 
Hilbert space can be mapped onto configurations of loops 
One remarkable feature of topological phases is that the 
ground-state wavefunction encodes many of the quasiparticle 
properties |2, 3, 4], which was exploited as far back as Laugh- 
lin's pioneering work on the fractional quantum Hall effect 
Therefore, the ground-state wavefunction plays a central 
role in the theory. Many of the properties of a QLG can be 
deduced by mapping the ground state to a classical statistical 
mechanical model. 

The toric code is the classic example of a QLG; the 
associated classical model is critical percolation. However, 
the toric code Hamiltonian is in a Z2 topological phase which 
is Abelian; i.e. all of its quasiparticles are AbeUan anyons. 
As we discuss below, non-Abelian topological ground states 
should be associated with critical 0(n) loop models with n > 1 
(The n = 1 case is equivalent to critical percolation.) 
However, local gapped Hamiltonians with these ground states 
are not known. In this paper, we show that such Hamiltonians 
do not exist. Hamiltonians with the desired ground states have 
been constructed in Refs. 9, i, but these models are gapless 
and describe critical points, not stable phases. It was conjec- 
tured that by perturbing such a critical model, one could drive 
the system into a gapped non-Abelian topological phase. In 
this paper, we analyze the instabilities of such critical models 
and show that perturbations fundamentally alter the nature of 
the ground state. For instance, one of the simplest relevant 



perturbations drives the system into the AbeUan I2 topologi- 
cal phase. Therefore, non-Abelian topological phases require 
more intricate Hamiltonians. 

Loop gas wavefunction.- QLGs can be realized in lat- 
tice models whose low-energy Hilbert space is spanned by 
states of the form \l), where the multi-loop L is an arbitrary 
collection of non-intersecting loops. Here we consider the 
'li-isotopy' wavefunction which generalizes the ground state 
(GS) of the toric code L7]: 
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In this expression, is the number of loops in the multi- 
loop L . The loops have a 'fugacity' d. In the toric code, d=\. 
This wavefunction has two key features which hint at its topo- 
logical nature: On non-trivial surfaces, there is a (degenerate) 
space of such wavefunctions corresponding, for instance, to 
different winding numbers. Second, the wavefunction ampli- 
tude is independent of the lengths of the loops. However, the 
latter is neither necessary nor sufficient for a gapped topolog- 
ical phase fs^. 

In a topological phase, the parameter d also determines the 
topological properties of excitations, as discussed in Ref. 0. 
Excitations can be studied by considering the ground state on 
a surface with punctures: each puncture can be viewed as a 
localized excitation which is specified by the boundary con- 
ditions at the puncture. If loops terminate at the boundary, 
the excitation is non-trivial (this is sufficient but not neces- 
sary). The amplitude to create a pair of such quasiparticles 
and annihilate them later is a measure of the number of states 
of such a quasiparticle called the quantum dimension. If there 
are quasiparticles with quantum dimension D, there will be 
^ degenerate states. For D > 1 and not an integer, there 
will be a large degeneracy which cannot be ascribed locally 
to the quasiparticles, so they will have non-Abelian braiding 
statistics. The universal properties of a topological phase are 
independent of any coordinate system; in particular, space and 
time can be interchanged. Therefore, the quantum dimension 
can be determined directly from the ground-state wavefunc- 
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tion. For a topological phase with ground state ([T]i, the quan- 
tum dimension of the fundamental quasiparticle is equivalent 
to the fugacity d. A loop can be viewed as the projection onto 
a fixed time slice of a pair creation and annihilation process. 
For d = V2, the fundamental quasiparticle has the same quan- 
tum dimension as the a-field in the Ising topological quantum 
field theory or the spin-1 /2 field in SU(2)2. For arbitrary fu- 
gacity d, the loop gas ground state ([T]) has no relation to any 
known topological phase. 

The fugacity d also determines the correlation functions of 
the associated classical statistical model, which is the 0(«) 
loop model with n = d^: 



biL) 
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Here, b{L) is the total length of the multi-loop L. For in- 
teger «, the right-hand-side is the expansion in powers of 
X of the Boltzmann weight e^P^ = ]!(,■)> (1 + xSi -Sj) for 
a model of classical interacting spins with 0(n) symmetry. 

Forx — n, Zo(n)(^) — (^cf'|1'o'')j so the equal-time ground- 
state correlations contained in the QLG's l^g'^) can be ob- 
tained from the known correlations of Zo(n) {x). For n<2 and 

x>Xc — n/ v2 + V5 — n, this model is in its low-temperature 



phase, which is critical (]9I]. It is necessary for the loops to be 
critical in order for jl'o'') to be the ground state of a topologi- 
cal phase; only then will the endpoints of a broken loop be de- 
confined. Therefore, for I'Pq'^) to be in a topological phase it 
is required that d < \/2- However, it is equally important that 
correlation functions of local operators decay exponentially in 
time since a topological phase requires a gap to excited states. 

Hamiltonian- The wavefunction ([TJ is the ground state 
of the following spin- 1/2 Hamiltonian, where the spins live 
on the edges of a honeycomb lattice: 
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The dominant first term (/ ^ K) enforces the constraint that 
the low energy Hilbert space is spanned by configurations 
with an even number of = 1 spins on edges e(v) around 
all vertices v. The loops are formed by the o~ = 1 edges. 
Fp = Yliep flips the six spins around a hexagonal plaque- 
tte p and P™ are projectors onto configurations with m loop 
segments around a given hexagon. (Here we define Fp so that 
it annihilates states with a single loop forming the plaque- 
tte boundary as in Fig. ^a); only configurations of the type 
shown in Fig. [Ttb) are not annihilated by P^^.) Notice that 
the Hamiltonian (|3]l includes processes shown in Fig. [na,b), 
but not those in Fig. [He, d). The most salient property of this 
Hamiltonian is that it is a sum of projectors which simulta- 
neously annihilate the wavefunction ^ which is, therefore. 
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FIG. 1: Four types of moves enacted by the plaquette flip term Fp. 
Term (a) enforces the loop fugacity d, (b) enforces isotopy invari- 
ance. Collectively they are known as 'd-isotopy' moves. Terms (c) 
and (d) are the surgery terms appropriate for d = I. 



the ground state. The low-energy spectrum of the Hamilto- 
nian (O is gapless iQ] for d < V2; the quantum dynamics 
is very inefficient at mixing states with different loop sizes, 
thereby resulting in gapless modes. 

For d — I a gap to all excited states can be opened with- 
out changing the GS \^q'^) by augmenting (O with the "loop 



surgery" terms in Fig. [nc,d) 
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This is the honeycomb lattice version of the toric code Hamil- 
tonian hence its GS is in a Z2 topological phase. Clearly, 
the existence of an energy gap above a ground state depends 
on the Hamiltonian. 

Augmenting the Hamiltonian ^ by the surgery terms of 
Fig. [nc,d) for d ^ I is not straightforward, as these terms 
generically do not conserve the number of loops and hence 
skew the loop amplitudes in Eq. ([T]). To preserve the correct 
amplitudes, one may propose a Hamiltonian 
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where A/ is the change in the number of loops when plaquette 
p is flipped. Although our Monte Carlo simulations show that 
this opens a gap A w 2K, the problem is that Al is a non-local 
operator: its eigenvalue depends on how the loop segments 
are connected away from a given plaquette p. 

Local Perturbations.- The question that we now want to 
address is: Can the non-Abelian state I'Po'') with d^lhe. the 
GS of a local, gapped Hamiltonian or is this specific to the 
Abelian case d =11 

One proposal is adding a local {k+ 1) -strand surgery term 
(not shown) given by a lattice version of the Jones-Wenzl pro- 
jector which also annihilates the GS \f d = 
2cosJt/(A; + 2). However, most of the states in this sequence 
d = 1, \/2, (1 +\/5) /2, . . . {k= l,2j3, . . .) occur at values of d 
for which the loops are not critical I9|]. The case d = \ {k=\) 
is the surgery term discussed above, and therefore we focus 
on li = V2 {k = 2). It can be argued that this term will not 
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FIG. 2: (color online) The local imaginary time correlation function 
C(t) for the loop gas with fugacity d and surgery term £ = 0.01. Data 
shown are for system size L = 16 and inverse temperature P = 16. 



open a gap, because there is a vanishing probabiHty for three 
long strands to meet [11]. 

However, it is possible for two strands to meet, so we con- 
sider the effect of the local 2-strand surgery term: 



^2 = ^0 



The state is no longer the GS of Hamiltonian (|6]), be- 

cause the second term is the surgery term for d — I. Nev- 
ertheless, one may hope that for e <C a weakly-perturbed, 
gapped GS will emerge which might still be in the desired 
non-Abelian topological phase. We simulate the model (|6]l on 
a torus of Lx L plaquettes with tilt angle 60° using a variant 
of the path-integral ground-state (PIGS) algorithm lfT2ll . We 
measure ground-state expectation values (0| A |0) by sampling 
the continuous-time path integral representation of 
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using local updates. We choose the projection time P large 
enough to project out all the states but the ground state. We 
find that this perturbation immediately opens a gap, since the 
time dependent local correlation functions 



C(T) 
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decrease exponentially as C(t) ^ A(T)exp(— tA), where the 
prefactor A(t) = (1 — exp{—zak^) / az is determined by a 
quadratic dispersion E{k) = A + ak^ above the gap A. The 
2-strand surgery term is expected to be a relevant perturba- 
tion, and an infinitesimal e should be sufficient to open a gap. 
As shown in Fig. |2] a substantial gap exists already for a small 
8 = 0.01. 
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FIG. 3: (color online) Length distribution of the longest loop in the 
Abelian (non-Abelian) quantum loop gas with loop fugacity d = \ 
(d — \fT). The multi-peak structure emerges from loops with differ- 
ent winding numbers. The data for various system sizes are rescaled 
by L''", with the Hausdorff dimension rf// = 7/4 {du = 3/2) in the 
Abelian (non-Abelian) case. 



We now need to determine whether this gapped ground state 
is in an Abelian or non-Abelian phase. While these two phases 
can, in principle, be distinguished by the fugacity of large 
loops, this is hard to measure. A simpler and more direct 
measurement differentiating the two phases is the Hausdorff 
dimension du of the longest loop. The Coulomb gas solution 
12] of the 0(n) loop model allows to calculate du exactly 1,1 ill 

with du = 7/4 for \^^^^) and dH = ^l1 for |^^^). We first 
measure these Hausdorff dimensions by sampling multi-loops 
il in a Monte Carlo simulation with weights given by the am- 
plitude d^^^^\ In Fig. [3] we plot the length distributions of 
the longest loops for d ^ \ and d = \/2. Rescaling the data 
for various system sizes by the expected Hausdorff dimen- 
sion we find an excellent data collapse. The characteristic 
multi-peak structure of the distributions originates from loops 
with different winding numbers j20l] . 

We have explicitly checked that the Hausdorff dimension is 
universal and constant for the full extent of a gapped phase by 
sampling the loop configurations created by domain walls of 
a ferromagnetic Ising model in its high temperature phase. At 
infinite temperature the Ising model on the triangular lattice is 
just percolation at the critical point, and the boundaries of per- 
colation clusters form a loop gas Although at any finite 
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FIG. 4: (color online) Data collapse of the rescaled loop length dis- 
tributions of various system sizes for small surgery term £ = 0.05. 
The loop fugacity isd = \fl, and inverse temperature P = 16. 



temperature above the critical point, coiTesponding to differ- 
ent GS wavefunctions within the same topological phase, the 
fugacity for small loops is changed, we find that the Hausdorff 
dimension of the longest loop stays li// = 7/4. 

Returning to model (|6]l, we now analyze the Hausdorff di- 
mension of the ground state in the gapped phase for small e. 
As shown in Fig. |4]the perturbation changes the Hausdorff 
dimension from c/// = 3/2 to t/// = 7/4 and the characteristic 
loop distribution to that of the Abelian phase. This can be un- 
derstood in a renormalization-group sense by noting that 
enforces the loop fugacity d only on microscopic loops, while 
the perturbation (which is the surgery term of the d = 1 loop 
gas) acts on large loops, and is thus a relevant perturbation, 
driving the system to the Abelian d — I fixed point. 

One might still hope that a different local perturbation to 
the Hamiltonian might open a gap but leave the system in the 
non- Abelian phase. We will now show that such hope is fu- 
tile by showing that they are power-law coiTelation between 
local operators in the GS \^\^^) for \ < d < Vl- While the 
correlations of most local operators are short-ranged, those of 
the plaquette flip operator Fp decay algebraically as r^^ with 
z « 3 for d — \/2 as shown in Fig. |5] and smaller values of 
z for d < A theorem by Hastings [ 14] then proves that 
because of this algebraic (and not exponential) decay of the 
coiTelation function between two local operators, the d ^ \ 
loop gas wavefunction cannot be the ground state of a gapped 
local Hamiltonian. 

The origin of the algebraic decay is the fractal nature of 
the loop gas: plaquette flips performing surgery between two 
segments of the same loop are correlated, since the change in 
loop number A/ is not just the sum of the changes of the in- 
dividual flips. Since the flip matrix element is d^' this results 
in an algebraically decaying congelation function for d ^ I. 
This argument applies not only to the wavefunction ([T]l but 
to any other wavefunction with long critical loops and a fu- 
gacity d ^ \ for large loops. Since the surgery term requires 
at least two loops to pass through the same plaquette or, in 
other words, four loop segments emanating from the same pla- 
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FIG. 5: (color online) Algebraic decay of the plaquette flip correla- 
tions between plaquettes at distance r in the QLG|1'')^'). 



quette, the {Fp{r)Fp(Q)) correlation function should not decay 
faster than the 0(n) model exponent associated with four-line 
vertices which gives z < 4 jg}]. 

We can thus prove that there exists no local, gapped 
Hamiltonian whose ground-state wavefunction is a loop gas 
with non-Abelian excitations. To construct local Hamilto- 
nians with non-Abelian excitations one needs to go beyond 
loop gases and consider more intricate Hamiltonians, such as 
string-net models (15, 3l. They can be constructed systemat- 
ically from loop gases with modified inner products lll7ll . 

For d — I the Hamiltonian Q is critical, but exhibits short- 
ranged equal-time correlations of all local operators. How- 
ever, by direct calculation we find that on-site operators have 
power-law decay in time, e.g., (o"'(0)O"'(t)) ~ 1/t. This is the 
firs t microscopic model exhibiting local quantum criticality 
iisll without dissipative baths. 
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lations were based on the ALPS libraries lUi 
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